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ABSTRACT 


This  paper  is  an  attempt  to  develop  a  method  for 
determining  the  slamming  pressures  and  loads  on  ship  hulls, 
when  the  relative  vertical  velocity  relationship  between 
the  wave  slope  and  keel  is  specified.  A  basic  two-dimensional 
theory  for  hydrodynamic  impact  is  utilized  which  avoids 
linearization  of  the  free  surface.  The  mathematical  model 
for  the  theoretical  developement  is  a  wedge.  A  method  of 
analysis  is  presented  which  yields  the  free  surface  shape 
and  pressures  associated  with  symmetrical  penetration  of  a 
fluid  surface  by  wedges  of  varying  deadrise  angle.  The 
anaylsis  is  based  on  the  principles  of  continuity  of  fluid 
flow  and  similarity.   The  principle  of  continuity  makep  use 
of  the  fact  that  the  volume  of  displaced  water  above  tfce 
original  water  line.   The  principle  of  similarity  accounts 
for  the  proportional  expansion  of  the  free  surface  with  con- 
tinued penetration  of  the  body.  The  free  surface  (consisting 
of  both  wave  and  spray)  is  considered  to  be  increasing  in 
such  a  manner  that  the  slope  and  velocity  of  the  surface  is 
determined  by  the  position  and  penetration  velocity  of  the 
impacting  body.  Once  the  correct  free  surface  shape  has 
been  determined  by  a  graphical  iterative  method,  for  any 
instant  of  time,  the  associated  velocity  potential  and 
pressure  can  be  determined  numerically  for  various  points 
on  the  body. 

Since  only  V-shaped  ship  sections  can  be  approximated  by 
wedges,  a  method  is  developed  by  which  the  pressures  on  a 
U-shaped  ship  section  of  a  low  deadrise  angle  can  be  calculated. 
Results  are  presented  which  resulted  from  calculations 
carried  out  on  a  ship  with  low  deadrise  angles  in  the  foward 
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sections*  For  comparison,  a  calculation  was  performed 
for  a  Liberty  ship  section  which  could  be  satisfactorily- 
approximated  by  a  straight  wedge  shape.  Also,  the  integrated 
pressure  load  was  determined  for  several  instants  of 
time  for  the  ship  being  considered. 

The  results  of  the  calculations  seemed  to  give 
realistic  values  for  the  pressure  distribution.  However, 
the  numerical  calculations,  as  developed  in  this  paper, 
are  quite  tedious.   Computer  solutions  need  to  be  developed. 
Also,  the  approximation  of  ship  sections  of  low  deadrise 
angles  by  blunt  nose  wedges  is  questionable,  since  the 
spray  root  thickness  has  to  be  approximated.   Therefore, 
the  free  surface  shape  in  the  spray  root  area  is  an 
approximation.   The  theoretical  results  should  be  closely 
checked  by  experiments. 
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I.    INTRODUCTION 

In  recent  years  there  has  been  a  distinct  trend  toward  increasing 
the  speed  and  displacement  of  ships.  Quite  often  a  ship's  speed  is  not 
governed  by  its  installed  horsepower,  hat  by  its  response  to  the  sea- 
way. Hie  term  'response    includes  not  only  the  ship's  motion,  but  al- 
so the  effect  on  the  ship's  structure  caused  by  the  motion.  One  possible 
result  of  a  ship's  motion  in  a  seaway  is  the  slamming  force  exerted  on 
the  forward  portion  of  the  hull  structure.  The  force  is  usually  caused 
by  the  emergence  of  the  bow  and  subsequent  impact  with  the  water  on 
reentry.  The  ship  may  vibrate  for  some  time  after  such  an  impact, 
and  plates  on  the  bottom  shell  may  be  damaged.  If  the  impacts  are 
severe  enough,  serious  structural  damage  may  occur,  and  even  if  the 
blows  are  of  small  magnitude,  but  frequent,  the  hull  will  possibly  be 
weakened  by  fatigue. 

The  earliest  theoretical  work  on  slamming  loads  resulted  from 
studies  on  the  impacts  sustained  by  seaplanes  when  they  landed.  The 
quantitative  evaluation  of  the  impact  forces  was  baaed  on  the    expand- 
ing plate  '  and    spray  root    theories.  The  expanding  plate  theory  ©rigi- 

e 
nated  with  Von  Karman  in  1923  (1)    .    By  this  theory,  the  water  flow 

pattern  around  a  wedge  penetrating  the  water  surface  is  taken  to  be  equi- 
valent at  each  instant  to  the  flow  about  a  flat  plate  of  the  same  width  as 
the  wetted  width  of  the  wedge.  From  this  analysis,  the  velocities,  accel- 
erations, and  resultant  pressures  and  hydrodynamic  reactions  are  ob- 
tained. Apparently  Von  Karman  *s  theory  is  good  for  low  deadrise  angles, 
but  is  only  a  crude  approximation  at  higher  angles  (2).  By  the  spray  root 
theory,  developed  by  Wagner  in  1322,  the  local  flow  phenomena  occuring 
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at  th«  edge  of  the  wetted  area  is  described.  The  high  pressures  that 
exist  near  the  edge  of  the  wetted  area  during  impact  are  connected 
with  the  formation  of  a  spray  jet  in  which  the  water  la  quickly  accel- 
erated to  a  high  velocity. 

.Seaplane  impact  theories  have  been  adapted  to  the  problems  of 
ship  slamming  (3,  4)  in  an  attempt  to  predict  the  pressure  response 
of  a  given  hull  form*  Experiments  have  been  devised  and  carried  out 
(4,  5,  6)  utilising  ship  models.  The  experimentations  have  been  carried 
out  under  conditions  allowing  three  dimensional  flow  about  the  ship  hull* 
which  are  obviously  the  conditions  met  by  ships  at  sea.   However,  three 
dimensional  flow  analysis  for  a  general  curved  bo$y  is  much  more 
complex  than  that  of  two  dimensional  flow.   Therefore,  theoretical 
work  has  been  primarily  concerned  with  the  flow  about  two  dimensional 
ship  sections  of  finite  length,  The  slamming  loads  for  the  individual 
sections  have  been  combined  by  strip  theory  to  obtain  the  pressure 
distribution  on  the  total  area  of  the  ship  bottom  affected  by  the  slam. 

There  Is  a  pressing  need  for  the  extension  of  the  available  two- 
dimensional  water  entry  theories.  In  the  expanding  plate  theory,  the 
free  surface  is  represented  by  a  aero  potential  line.  The  impact  of 
ship  hull  shape  bodies  on  a  water  surface  will  usually  result  in  the 
formation  of  a  curved  free  surface,  and  a  spray  root.  Theoretically 
the  presence  of  a  highly  curved  free  surface  may  not  be  represented 
by  a  straight  sero  potential  line.  It  seems  likely  that  the  tv9  surface 
is  nearly  a  horisontal  line  only  for  bodies  of  high  deadrise  angles*  and 
mat  the  spray  root  generated  at  low  deadrise  angles  disputes  the  assump- 
tion of  a  horisontal  aero  potential  line.  In  general*  for  any  of  the    fit- 
ting   techniques  ( ellipse  and  circle  fitting  as  well  as  the  expanding 
plate),  the  agreement  between  theory  and  experiment  has  been  rather 
poor  (2)*  even  for  simple  wedge  shaped  bodies.  Therefore,  a  more 
exact  solution  ot  the  hydrodynamic  impact  problem  would  be  desirable. 
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The  object  of  this  work  is  to  formulate  a  basic  two-dimensional 
theory  for  hydroctynamiC'  impact  which  avoids  linearisation  of  the  free- 
surface  boundary  condition,  and  which  can  be  adapted  to  the  ship  im- 
pact problem.  The  theory  is  baaed  largely  on  the  work  of  Fierson  (7) 
who  developed  an  iterative  numerical  solution  for  a  wedge  entering 
the  water  at  constant  velocity.  While  a  wedge  is  used  as  the  mathe- 
matical model  in  this  paper,  the  velocity  is  not  considered  to  be  con- 
stant. 

The  problem  of  determing  slamming  pressures  actually  consists 
of  two  parts;  first,  the  dee  cription  ot  the  ship  motion  in  a  seaway  which 
leads  to  slamming,  and  second,  the  description  of  the  pressures  and 
forces  on  the  hull  caused  by  slamming.  Ship  motions  are  not  treated 
here,  but  an  expression  for  ship  motion  is  needed  tie  determine,  at  any 
instant,  the  relative  vertical  velocity  of  the  ship  with  respect  to  the 
water  surface.  Such  an  expression  'was  obtained  from  experimental 
work  for  a  Liberty  type  ship,  discussed  in  reference  (3). 

Because  of  the  difficulty  of  analytically  describing  the  flow  about 
an  arbitrary  convex  bo<fy,  the  expressions  for  the  velocity  potential  de- 
rived in  the  paper  are  based  on  the  water  entry  of  a  wedge.  The  dead- 
rise  angle  of  the  wedge  is  taken  to  be  approximately  the  same  as  the 
deadrise  angle  of  the  ship  section  under  consideration,  While  the  above 
assumption  appears  to  be  somewhat  gross,  in  the  ligfci  of  past  experi- 
mental and  theoretical  work  it  does  not  seem  to  he  bad.  Both  expert 
mental  and  theoretical  calculations  have  shown  that  the  highest  trans- 
verse pressure  distributions  occur  on  a  given  ship  section  shortly  after 
it  hits  the  water  (3,  @,  8).  The  theoretical  calculations  were  based  on 
the  expanding  plate  theory,  with  a  modification  for  the  hull  section  cur- 
vature and  added  mass  of  water.  However,  the  portion  of  the  hull  that 
experiences  the  largest  pressure  loading  is  the  bottom  portion  which, 
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for  most  ships*  is  roughly  ©quivaiezit  to  a  wedge  of  low  deadriae  angle, 
or  can  at  least  b«  approximated  {at  leaut  for  the  period  just  oaf  or*  and 
after  slamming)  by  a  wedge  of  a  specified  angle.  It  is  not  sufficient 
just  to  compute  me  pressure  distribution  at  ih@  instant  of  slamming, 
but  pressure  distributions  for  several  instants  before  and  after  slam- 
ming are  also  require©,  for  a  complete  picture  of  me  pressure  loading. 
From  a  structural  viewpoint,  it  is  the  pressure  loading  and  its  time 
variation  *aich  are  important,  and  not  peak  pressures  which  last  only 
for  a  few  micro -seconds. 

The  development  of  the  theory  is  presented  in  Chapter  DL  The  re- 
sults of  the  theory,  as  applied  to  a  given  ship  for  which  other  theoretical 
and  experimental  data  are  available,  are  presented  and  discussed  in 
Chapters  IB  and  IV. 
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II,    PHGCEDU&B 

A.  General  The  theoretical  calculations  are  carried  cut  for  a  ship 
approxlsoaiely  the  sis*  of  a  Liberty  ship.  The  ac&ial  ship  li&es  were 
calculated  from  the  scries  ;s  for  which  Protease  ,  If.  Lewis 
obtained  his  values  of  inertia  coefficients  <S)„  Protestor  J.  B.  Kerwin 
obtained  the  offsets  and  polynomial  expressions  for  the  ship  lines  by 
the  use  of  the  ?®$Q  IBM  computer.  Also.,  a  calculation  is  performed 
for  comparison  for  a  section  of  the  Liberty  ship  nasd  in  references 

(3,  4). 

B.  Determination  of  &htp  Motion      The  coordinate  sjmimss.  xatod  for  the 
determination  of  ship  motion  is  the  same  one  used  la  reference  (10), 
The  eoor<&aaie  origin  is  taken  to  fee  at  midships,   I  i  fc  y9  s)  coordi- 
nates represent  linear  displacements  of  the  center  of  gravity  along 
fixed  axes.  In  this  paper,  only  heaving  and  pitching  motions  are  con- 
sidered* therefore  only  the  a  coordinate  is  of  interest.  The  s  axis  is 
vertical  and  s  is  positive  upward.  Of  the  three  angular  displacements 
(roll,  pitch,  and  yaw),  only  the  pitch  angle  is  of  interest.  Positive 
indicates  bow  up.  Also,  the  <X9  Y*  Z)  coordinate  system,  which  is  re* 
lated  to  the  ship's  principal  axes,  is  introduced.  The  Z -coordinate  is 
positive  up,  the  X -coordinate  positive  forward,  asd  the  Y  -coordinate 
positive  to  port.  The  pertinent  coordinates  and  their  positive  directions 
are  shown  in  Fig.  1. 

In  order  to  obtain  a  complete  picture  of  the  afcip  m  at  any 

instant  of  tima,  it  is  accessary  to  know  the  instantaneous  draft  and  re- 
lative velocity  between  the  keel  and  water  surface  for  any  station. 

The  pitching  motion  is  assumed  to  be 

y   *  y       cos  (die*  ♦  €  ) 
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FIGURE  I    PRINCIPAL  COORDINATES  AND  AXES  OP 
REFERENCE 


FIGURE  II   RELATIVE  ELEVATION  OF  KEEL  ABOVE 
DISTURBED  WATER  SURFACE 
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and  the  heave  is  assumed  to  be 

s  «  s     cos  (fi^z-t  +  6  K 
us 

The  wars  elevation  at  any  station  along  the  ship  is  given  by 

'  m  a 

where  oj<.  •         <v  +  -£*  )   *  period  of  encounter  between  ship  and  waves 
in  head  seas. 

The  relative  position  of  the  keel  above  the  disturbed  water  surface 
is  given  by: 

»jt-r-»-XH'  tH  U) 

for  any  station  (see  Figure  II).  Emergence  is  associated  with  *    <  0, 

The  relative  velocity  between  the  ship  and  wave  at  a  distance  X 
from  amidships  is  given  by 

ix  »  r  -  i  -  x  f  (2) 

and  CU0  is  the  absolute  angular  frequency  of  the  wave  and  equals   -~£- , 

1 1  m 

andu;*    ~x    •  -^     ccfe  aia(fcfct  +  €). 

1        el  '  m 

From  reference  (11),  experimental  values  of  3  and     were  obtained 
for  a  5.  ft  ft.  model  of  a  Liberty  ship  in  regular  waves  of  a  specified 
height.  The  values  given  for  the  model  are  the  same  for  the  full  scale 
ship,  with  the  exception  of  the  heave  amplitude,  which  is  sealed  up  for 
the  ship  by  a  ratio  of  the  lengths.  The  pertinent  variables  for  the  mo- 
del, full  scale  ship,  and  equivalent  Lewie  ship  are  tabulated  in  Table  I. 

The  calculations  of  the  relative  velocity  between  the  keel  and  the 
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water  by  Bq„  (2)  will  result  in  &  sinusoidal  variation  of  velocity  for 
any  station.  While  this  is  not  correct*  it  is  shown  in  reference  (3) 
that  this  is  a  good  approximation  for  the  relative  velocity  of  a  ship 
subject  to  slamming  motion  in  regular  waves.  In  other  words,,  the 
change  in  relative  velocity  due  to  slamming  loads  is  very  small* 

Slamming,  as  generally  defined*  is  the  sudden  change  of  the 
acceleration  of  the  ship.  While  the  relative  velocity  does  not  deviate 
much  from  that  predicted  by  (2),  the  influence  of  slamming  on  the  accel- 
eration is  quite  marked  (3).  However,  in  this  study  of  slamming,  the 
pressures  are  governed  by  the  velocity  and  penetration  of  the  keel  in- 
stead of  the  acceleration.  For  the  ship  hulls  under  consideration,  the 
portion    of  the  ship  between  stations  9  and  a  Is  subject  to  slamming. 
From  reference  (3)  the  instant  of  slamming  is  determined  to  be  in  the 
vicinity  of  t  *  1. 29  sec.  The  penetration,  and  relative  velocities  of  sta- 
tions 0  through  5  are  graphed  for  the  interval  from  t  •  1. 2$  sec.  to 
t  »  1,44  sec.  in  Figs.  HI  and  IV.  These  values  will  be  used  to  deter- 
mine the  slamming  pressures  on  the  forward  portions  of  the  ship.  The 
values  used  to  plot  the  curves  in  Figures  IH  and  IV  were  computed  from 
equations  ( la)  and  ( lb)  which  resulted  when  the  amplitudes  and  phase 
angles  from  Table  I  were  substituted  into  equations  (1)  and  (2). 

»    ■  -12.  $  cos  {coet  +  H| )  *  g» 3$  cos  (ftfct  -0.  ©23)-0. 122Xcos<<^t+0. 26 

(la) 

■^  •  *9. 3  sinfcvt*  !fr  >"s- u  8in  ****  *"  •  •*•>-  •  i3*  ^^ ***  *  •  ***> 

<lb) 

C.  Theoretical  Analysis 

Conditions  of  the  Problem 
1.        The  problem  is  treated  as  a  two-dimensional  one.  The  velocity  po- 
tential and  pressure  distribution  will  be  calculated  for  two-dimensional 
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TABLE    I 


Model 

Liberty 

ship 

Lewis 

type  ship 

Length  (BP)      L,   ft. 

5.5 

416/0 

420.0 

Beam,    B,   ft. 

0.  732 

56.9 

57.5 

Draft,   H,   ft. 

0.21 

15.9 

16.0 

Scale  factor 

1.0 

75.7 

76.5 

Speed  V,   knots 

1.1 

9.6 

i 

9.7 

Wave  length       X  ,   ft. 

5.0 

379.0 

382.0 

Wave  height        h,   ft. 

0.334 

25.3 

25.5 

Pitch  amplitude     lf>      ,   degrees 

i  m 

7 

7 

7 

Heave  amplitude     Z     ,   ft. 

0.007 

5.30 

5.35 

Pitch  phase        £    ,    degrees 

14.44 

14.44 

14.44 

Heave  phase         & ,   degrees 

-35.7 

-35.7 

-35.7 

■r-,                       »                       ,  .      rad. 
Frequency  of  encounter  cU?,  

8.74 

1.01 

1.00 

Period  of  encounter     T   ,    sec. 

e 

0.72 

6.28 

6.28 

rad. 
Wave  frequency       QJ^,    

6.36 

0.73 

0.73 
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sections  (each  *  station  spacing  in  length),  and  then  combined  by  strip 
theory  to  yield  the  pressure  load  on  the  forefoot  of  the  ship. 

2)  Each  impacting  section  is  assumed  to  he  perfectly  rigid  and  to  be 
symmetrical  about  the  vertical  c  eater  line  plane. 

3)  The  velocity  of  penetration  is  not  constant,,  but  varies  according 
to  equation  (2). 

4)  The  flow  field  is  assumed  to  be  irrotational  (a  necessary  condition 
for  the  existence  of  a  velocity  potential). 

5)  The  fluid  is  assumed  to  be  frictioaless  and  incompressible,  and 
the  velocity  of  penetration  is  assumed  to  be  high  enough  for  the  effects 
of  gravity  and  surface  tension  to  be  neglected. 

6)  On  the  basis  of  the  assumption  of  irrotational  flow,  the  particles 
of  the  original  surface  remain  on  the  surface  and  no  new  particles  are 
added  to  the  surface  (12).  Since  the  fluid  is  incompressible,,  the  arc 
length  along  th«  surface  from  the  impacting  body  to  any  given  particle 
is  constant  and  independent  of  time. 

The  preceding  conditions  do  not  define  the  shape  of  the  free  sur- 
face. However,  for  a  wedge  shaped  body,  which  is  the  mathematical 
model  used  for  the  theory,  two  general  conditions  can  be  formulated 
which  will  define  the  shape  of  the  ttm  surface; 

*)  Continuity;  The  fluid  displaced  by  the  body  must  appear  above 
the  original  water  surface  in  the  form  of  a  wave  and/ or  spray.  The  arc 
length  of  the  surface  must  be  constant. 

W  Sinat  larlty s     Since  the  impacting  body  is  symmetrical,  and  the 
immersion  is  normal  to  the  surface,  the  field  of  flow  is  divided  into  two 
halves  by  the  centerline  of  symmetry.  By  neglecting  the  gravity  forces, 
the  surface  shape  and  tynamlc  state  of  the  field  of  flow  can  be  defined 
entirely  in  terms  of  the  penetration  and  velocity  of  the  body  (?).  If  the 
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penetration  velocity  is  constant,  the  flow  field  and  surface  shape  must 
be  geometrically  and  dynamic  ally  similar  at  all  times  a  However,  in 
the  present  case,,  the  penetration  velocity  is  not  constant,  and  there- 
fore the  surface  shape  is  assumed  to  be  only  geometrically  similar  at 
different  times. 

Application  of  General  Conditions 

CONTINUITY      The  general  shape  of  the  free  surface  due  to  penetra- 
tion by  a  wedge  shaped  bo4y  ia  shown  in  Fig.  V.    As  shown  in  the  figure* 
the  free  surface  can  be  thought  of  as  being  made  up  of  spray  and  a  wave. 
The  two  merge  together  in  the  "spray  root'  area.  Ia  determining  the  free 
surface  shape,  the  wave  rise  {$)  ia  given  by  equations  11  and  13  of  refer- 
ence (13). 

.  LSJ85JL  ,  t  „»-»    c_   .   j  ( 

-  c  y 

This  is  the  same  expression  which  gives  the  surface  rise  in  the  expand* 
ing  plate"  theory.  The  watiad  half  width  (c>  is  the  point  on  the  body  which 
is  wetted  by  the  piled  up  water  in  the  absence  of  spray.  The  wetted  half 

width  at  the  original  water  level  is (Wagner;  see  reference  14). 

The  continuity  condition  is  satisfied  If  Area  A  *  Area  B,  and  if  s  a  y. 
However,  the  continuity  condition  Itself  is  not  sufficient  to  determine 
the  free  surface  shape  since  any  thickness  of  spray  may  be  associated 
with  an  unlimited  variety  of  waves,  and  still  not  violate  the  continuity 
requirement.  Also,  the  free  surface  in  the  spray  root  area,  where  the 
surface  wave  and  spray  merge  together,  could  assume  any  one  of  many 
shapes. 

Similarity    The  principle  ©f  similarity  for  the  proportional  increase  of 
the  wave  and  spray  with  continued  penetration  is  the  second  considera- 
tion. Not  only  must  the  displaced  fluid  appear  above  the  original  water 
surface*  but  the  wave  and  spray  must  be  increasing  in  such  a  manner 
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that  tht  shape  and  velocities  Mm  everywhere  determined  fey  the  position 
and  penetration  velocity  of  the-  entering  baity,  By  a  geometrical  con- 
struction and  integration  process  proposed  by  Wagner  (see  reference  7), 
it  ie  possible  to  relate  the  surface  shape  and  velocities,  so  that  for  any 
spray  thicknesa,  one  corresponding  wave  and  spray-root  shape  can  be 
determined  which  meets  the  requirements  of  continuity  and  similarity. 
However,  Wagner's  proposal  was  tor  a  constant  entry  velocity,  the  follow* 
ing  theory  is  developed  for  1  variable  penetration  velociv 

In  an  ideal  fluid  where  gravity,  viscosity,  and  surface  tension  are 

neglected,  the  snap*?  if  the  free  surface  is  uniquely  determined  by  the 

shape  of  the  impacting  body,  and  the  sise  of  the  shapes  is  determined 

e 
by  the  amount  of  penetration  . 


Consider  the  several  stages  of  the  penetration  process  of  a 
shaped  bo^jr  im  Figure  VI,  using  an  y~«  coordinate  system.  The  penetra- 
tion of  the  body  into  the  fluid  eaaflas  the  initially  flat  surface  to  rise  as 

a  wave*  and  the  nearest  particles*  are  deflected  out  along  the  wedge  sur- 
face as  spray,  The  resulting  surface  shape  expands  as  the  wedge  pane* 
trates  deeper.  The  resulting  surface  shapes  are  o*j  aAbt,  o^s*  tfci  etc. 
For  the  same  relative  position  ia  the  fluid  and  on  the  surface,  the  prm~ 
•-are  and  velocity  depend  upon  the  surface  ahape,  and  must  bear  the  same 
relationship  to  the  penetration  velocity  at  any  instant  of  time  <?)♦  The 
above  statement  applies  to  the  relative  position  in  the  fluid,  end  ®®*  to 
a  particular  fluid  particle.  The  »  m  are  indicated  by  o\  ,  »*« ,  I 

and  a*  *  a*  ,  **,  ate.  If  the  penetration  velocity  was  constant  for  the  body, 
the  velocity  of  the  fluid  at  corresponding  points  would  be  equal.  That  is 

Of  course  this  statement  is  only  true  for  a  certain  definite  time  inter- 


val after  impact.  As  ^he  bo%  penetrates  deeper,  much  of  the  free  sur- 
face shape  will  be  dissipated  as  spray,  and  as  the  penetration  i 
decreases,  the  influence/  of  gravity  will  be  felt. 
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U    »  U     •  U     »  K    *      where  K    represents  a  proportionality  factor 

e\      a.jj       a«?       ax  a 

between  the  velocity  U  at  positions  %  ,  •*,  atc.p  and  the  penetration 

velocity,  a  .  However,  for  varying  velocity,  U     j  U     ^  U     ,  but  the 

*  »i        *«?        a$ 

velocity  at  corresponding  points  wiU  still  be  proportional  to  the  pene- 
trating velocity. 

VS^'  °b>  •*!,*,<«•>  (4> 

However,  it  la  not  unreasonable  to  assume  that  K     «  K      a  K     ,  K.   * 

84        a*        a#       oj 

K    »  K.    ,  etc.    Therefore,  the  velocity  at  corresponding  points  will  be 
directly  proportional  to  the  instantaneous  penetrating  velocity. 

The  velocity  at  any  point  on  the  free  surface,  such  as  a?  or  b«, 
is  determined  from  consideration  of  this  condition  of  similarity.  Con- 
sider point  h~:  for  example,  A  particle  on  the  surface  at  point  b©  is 
instantaneously  moving  along  the  radius  with  a  velocity  U  .  and  at  the 

|SJ 

same  time  sliding  past  the  radius  toward  the  spray  region  with  a  velocity 
U       (see  Fig.  VII). 

The  vector  sum  of  the  two  motions  ( U      +  0^  )  yields  the  instan* 

9%  Sfe 

taneous  resultant  velocity  U    of  the  point  in  question. 

The  radial  expansion  velocity,  U  ,  the  velocity  of  the  surface 
past  the  radius,  U  ,  and  the  resultant  velocity,  U,  are  each  directly 

proportional  to  the  penetration  velocity,  z  .  Therefore,  for  constant 

x 

entry  velocity,  or  for  an  average  velocity  taken  over  a  small  time  in- 
terval, the  velocity  components  ( U  ,  U  )  at  a  position  are  constant,  and 

r      s 

the  vector  sum  may  be  obtained  from  the  radial  distance  and  the  arc 
length. 

Consider  the  point  b* ,  on  the  curve  in  Figure  VI.  The  radial 
distance  r  is  given  by: 
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The  vector  9  is  constructed  tangent  to  the  surface  at  br.   so  that 
s  equate  the  arc  length  $'h^  ^ 

"1   »      /  U"     dt 
The  vector  earn  ia  then: 

-  -   c 

s   +   r   *  U.  dt  (5) 

The  individual  fluid  particles  do  not  follow  the  radiating  lines, 
but  remain  at  a  fixed  arc  length  from  0'  .  The  motion  of  a  particle  ie 
therefore  a  function  of  its  relative  position  on  the  free  surface,,  and  is 
given  by  the  following  integral  relations 

S   -h        f     •  <«a) 

Jo 

A 


•  .-.  Woi 


>: 


Equations  (0)  are  referred  to  the  y  and  a  axes  in  Figure  VI.  The  motion 
in  to.  y  *r.ctio»  In  •  Dm  tofrvl  t  to  «pr«M«  by  ?    -  £  .  £  M« 
the  original  position  of  the  particle.  The  symbol  3   represents  the  vert- 
ical elevation  of  a  particle  (a -direction).  Now  the  velocity  of  a  particle 
(v,  w)  is  to  be  determined  for  any  point  in  its  path,  so  that  Equation* 
(6)  can  be  solved.  There  does  not  seem  to  be  a  simple  analytic  solution, 
but  a  graphical  solution  can  easily  be  carried  out.  However,  the  solution 
can  be  simplified  by  changing  the  variable  of  integration  from  time  to  are 
length.  This  is  done  in  the  following  manner.  If  the  penetration  velocity 
was  constant,  the  arc  length  and  times  at  two  corresponding  points  on 
the  free  surface  could  be  related  fays 

4GU  tt  til 

M  ti  ti 

(reference  (?) ).  However,  for  varying  velocity  at  time  t,  V  .  •  K  .  1  ft*  >» 
where  K  fe  is  the  proportionality  factor  which  relates  ihe  velocity  of  the 
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surface  sliding  past  the  radiu  penetration  velocity,  How  *n 

the  time  Interval  ft  mx  assume  that  U  .is  constant  and 

Mi 

equal  to  its  average  value.  According  to  reference  ( 15),  an  average  valve 
for  a  function,  ftt),  can  be  defined  as  follows: 


^f(t) 


t  (t+  ^  )  ♦  f  (t  -  ^~-> 


tt    ♦  t*> 

Therefore,  at  time  -3a— a. ,  th«  average  value  of  V  .  can  be  assumed 

to  be  constant  over  the  interval  t%  to  t© .  For  any  Uaia,  t-j  , 


« 


a 


(k    1  '        t 

By  this  reasoning,  an  expression  for  arc  length  in  terms  of  time  can 
be  obtained  for  any  time  interval,  At  : 

t.__ 

where  t  *  the  earlier  time,  and  S  »  corresponding  arc  length  on  the  sur- 
face at  time  ta. 

Mi 

III  \  dt  « ■   -  -gy-     da 

Substitution  for  dt  in  Equation®  <6a)  and  (6b)  yields 

$     -I'    fvC"-^1-  d»>  "  S.  (%-  d*  <?a) 

/-OO  ■*> 

I  •       S0\  ^  dm  ilh> 

Jo 

Equations  (7a)  and  (7b)  can  be  solved  numerically  as  follows: 


Referring  to  Fig,  VI,  it  is  seen  that  the  penetration  velocity  has 
approximately  a  linear  variation  over  the  time  interval  under  considera- 
tion. Since  the  particle  velocities  at  corresponding  points  are  directly 
proportional  to  the  penetration  velocity,  they  can  also  be  assumed  to 
vary  linearly.  Using  this  assumption,  the  integral  in  Equation  (5)  can. 
be  expressed  as 
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"S  +  r  *  J  Ub  dt  «  Ufet  (5) 

where  U.   is  the  average  value  of  U.   for  the  time  interval  under  question. 
From  this,  the  average  value  of  the  velocity  components,  in  Equations 
(7a)  and  (7b)  can  be  found  (vt  and  wt).  Therefore,  a  step  by  step  inte- 
gration of  Equations  (7a)  and  (7b)  can  be  carried  out  as  shown  in  Appen- 
dix A.  Any  trial  free  surface  shape  which  satisfies  continuity  can  be 
used  for  the  determination  of  the  surface  velocities.  The  velocities 
are  then  used  to  determine  a  new  surface  shape.  This  procedure  is 
carried  out  until  there  is  a  negligible  difference  between  assumed  and 
computed  values  of  surface  shape. 

The  Integral  equations  (7a)  and  (7b)  relate  the  surface  rise  at  a 
point  to  all  the  points  more  distant  from  the  body.  Therefore,  the  path 
of  each  particle  an  the  surface  follows  the  same  related  path,  and  the 
conditions  for  similarity  are  satisfied.  It  would  appear  that  the  integral 
equations  reduce  the  number  of  spray -wave  combinations  to  one  for  each 
spray  thickness  (7). 

The  spray  root  thickness  itself  was  shown  in  reference  (7)  to  be 
related  to  wetted  width  (c)  and  deadrise  angle  (3): 

6  •  c  (  4-  )»  I  \ 

Potential  on  Free  Surface 

Referring  to  Figure  VII,  the  tangential  velocity,   ~x   along  the 
surface  at  a  position,  b  is  given  by  -  U  .  +  U  .   cos  Y  .  For  points  along 
the  surface  (at  a  specific  depth)  eo®T=     -r-    *  . 


Therefore, 


[  •    B  u  -  U*.  cosY  *  U  -    -  U  K   ~  (9) 

Ss         tli        rfo  rb    ds 


+)   dr  dr  , 

~  for  a  particular  surface  should  not  be  confused  with  -       *    /s  for 

a  constant  position,  such  as  b. 
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Since  Ufife  and  0^  are  proportional  to  i  „  and  *    can  be  assumed  equal 

to  a  mean  value  over  the  time  interval  under  consideration,  U    and  U 

e  r 

can  be  considered  approximately  constant  over  the  eame  time  interval. 

A  differential  equation  can  be  obtained  by  multiplying  Sanation  (8)  by 

tda.  (tUe^s,    tU  «r): 

r 

~  td<£  •  sds  -  rdr 

Integration  yields 

IsF  -  v*"  I 


which  is  an  approximate  expression  for  the  velocity  potential  on  the  free 
surface  with  varying  penetration  velocity., 

Potential  and  Velocity  Distributions  on  Body  Surface 

Assuming  that  the  free  surface  shape  for  a  given  deadrise  angle 
is  unique,  end  is  a  result  of  potential  flow,  the  velocity  and  potential  at 
any  point  can  be  represented  by  a  complex  analytic  functions 

The  Cauchy  integral  equation  may  be  applied  to  the  complex  func- 
tion as  described  in  reference  (Th  By  Cauchy' s  integral  formula,  a  value 
of  an  analytic  function  in  a  region  t&  determined  by  its  values  on  the  boun- 
dary. There  is  no  choice  of  ways  in  which  the  function  can  be  defined 
at  points  away  from  the  boundary  once  the  function  is  defined  on  the  boun- 
dary. Every  change  of  values  of  the  function  at  interior  points  must  be  ac- 
companied by  a  change  of  its  values  on  the  boundary  if  the  function  is  to 
remain  analytic  ( IS).  The  integral  formula  for  the  point  within  the  boun- 
dary is; 

"Via   /^r-  <«  «» 
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s  is  the  complex  coordinate  *f  the  position  in  the  field.  Since  poiate  on 
the  solid  boundary  are  being  investigated,  the  bounding  contour  circles 
only  halfway  around  the  point  so  thai  Set  becomes  si.  Also*  #s)  ■  a  +  iv 
*  UU)  on  the  boundary*  and  equation  (It)  becomes: 

The  expression  can  be  separated  into  the  real  *&&  imaginary  parts  by 
placing  the  origin  of  coordinates  at  the  point  s    under  consideration, 
and  writing  a  in  the  polar  forms 

Substitution  of  the  expression  for  %  and  ds  into  Equation  (1$)  yields 

o  si  J  l  J 

Rearranging  terms  give* 

tru  +  av-    jd  (laSi>  *  v*te   +  i  £ud<*    -vd(laR)"] 

so  that  tie  real  and  imaginary  parts  may  be  separated; 

v    «   i  /  vd6   +  0d(inE)  Ufa) 

o         > 

U    *i    I    d^    -vd(laR)  (1Tb) 

o      *  J 

The  construction  for  performing  the  numerical  integration  of  the 
above  equations  is  shown  in  Figure  Yin*  The  a*  linates  were 

oriented  so  that  the  y~exia  (Imaginary  axia)  is  along  the  wedge,  positive 
toward  the  spray  tip.  The  x-sxis  {real  axis)  if-  normal  to  the  body  and 
positive  out*  The  potential  on  the  bod(r  surface  may  be  obtained  by  Inte- 
grating the  velocity  along  the  wedge  from  the  spray  tip  to  the  apex, 

Pressure  Distribution 

The  general  equation  for  fluid  p&%»mr«t  in  a  potential  flow  of  in- 
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compressible  fluid,  is  used  to  determine  the  pressure  distribution 
on  the  body  (13) 

£  •  -  r(t>  (18) 

The  gravity  fore*  (Q)  and  extraneous  Impulses  <F  (t) )  ere  aero  sine* 
gravity  haa  bean  neglected,  and  the  fluid  was  assumed  to  be  at  reat  at 
Infinity.  Therefor*  Equation  (  18)  becomes 


a 


p  V 

r     IT     I  *    ' 

The  resultant  field  velocity  (U)  can  be  determined  on  the  free  surface 
(Equation  ($)  )„  and  on  the  body  (Equation  f  11)  L  Now  it  remains  to  de- 
termine an  expression  for 

Assume  that  the  depth  of  penetration  at  time  t%  is  given  by  the  aver- 
age penetration  velocity  <«  >,  in  the  interval  t  *  o  to  t  *  t4 ,  multiplied 
by  t4  .  This  la  the  same  aa  assuming  constant  velocity  penetration  in 
the  specified  time  interval*  At  any  position,  (y*  f  %J     in  the  fluid*  the 
potential  is  therefore  a  linear  function  of  time  in  the  interval  to  %%  (7)« 

For  any  time  interval  %  to  %  assuming  i. 

However,  for  a  fixed  point  (y,.  s)  in  the  flow  field,  the  above  expression 
i«  not  correct.  If  the  point  <y,  »)  coincides  with  the  position  r t  ,  a      at 
tx ,  at  time  !•%#!  t  the  poatUon  r%  +4t)J  8    is  a  distance  &r  * 
from  (yf  a).  For  smell  At,  the  resultant  vector  velocity  can  be  assumed 
constant  tw  the  region. 

At  t  *    .    4>(y,  a)  *  4>r        and  at  t  *  t4  ♦  *-t  4*y#  *>t  +At    *^t  "  Uco**(   r 

See  Figure  IX,  (T    a  )  refers  to  the  distance  to  the  same  relative  point 
in  the  fluid  at  a  deptfrof0  peneiratioa  a  .  Therefore,  as  time  Increases,  T% 
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FIGURE  IX   FIXED  AND  RELATIVE  POSITION  POINTS 
IN  THE  FLUID 
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U  cos  a  u»  the  velocity  component  elong     .  Therefore, 

■/•*>t,  f4t        (*.*>&*<  rcoea      r      (21) 


From  Equation  (2&) 


■ 


t»    ♦  &t 

V»        *     A  _  ..*«—__       IK        ft)  A      . 


Station  (21)  then  becomes 

A  * 

A*  »  •  vv     •*  tt&t  cos  a  (22) 

Taking  the  limit  of  Equation  (22)  as  At  approaches  0: 

n 

(23) 

Equation  (22)  now  applies  to  ail  points  in  the  fluid.  The  pressure  distribution 
if  given  fey  substituting  Equation  (23)  in  Equation  (1$). 

JL   «    .  ^.  4.   ^      u  cos  *  *  «-.  (24) 

On  the  body  surface,,  the  pressure  may  be  found  by  using  the  values  of     , 
-4.  U,  and  a  already  determined.  The  total  load  on  the  body  is  obtained 
from  the  integration  of  the  pressure  across  the  projected  wetted  width. 


Application  to  Ship  Sections 

There  are  twe  possible  methods  of  solving  the  water  entry  problem 
of  bodies  of  general  shape  by  considering  the  free  surface  shape.  The 
first  is  to  construct  the  free  surface  based  on  the  boundary  ot  the  actual 
ship  section  as  it  enters  the  water.  However,  this  method  is  questionable 
since  the  similarity  principle  has  not  been  proved  to  hold  for  a  body  of 
general  shape.  The  second  method  is  to  approximate  the  ship  section  by 
a  wedge-like  shape  of  an  equivalent  de&drise  angle. 

The  period  of  interest  in  a  slamming  investigation,  for  a  certain 
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•hip  section,  is  the  time  period  immediately  after  impact  with  the 
water,  The  maximum  impact  force  is  produced  either  at  the  instant 
of  Impact,  or  immediately  after  entry .  For  the  Liberty  ship  ixvmti- 
gated  in  reference  (3),  It  was  found  that  the  slamming  pressures  were 
of  little  significance  by  the  time  the  draft  at  station  &  was  9,2ft.  In 
fact,  the  time  interval  ot  eigaif leant  pressures  was  l»m  than  0,1?  sec. 
For  this  problem,  the  time  interval  w&to®r  investigation  is  0. 24  sec. 

W&r  a  hull  section  with  sharp  curvature,  the  spray  will  act  neces- 
sarily adhere  to  the  boundary.  From  photographs  of  wedge  impacts  in 
references  (14,  17),  the  spray,  tiler  it  leaves  the  wedge  surface,  shoots 
off  at  an  angle  roughly  equivalent  to  the  deadrise  angle  (see  Fig.  X). 
Therefore,  it  99emu  possible  that  v- shaped  sections,  or  sections  with 
fine  lines,  could  be  approximated  by  a  wedge  of  a  certain  deadrise  angle. 
However,  slamming  is  more  dangerous  for  sections  with  small  deadrise 
angles,  It  is  obvious  that  these  sections  could  not  satisfactorily  be  re- 
presented by  a  V-  shaped  wedge. 

For  U-shaped  sections  with  small  deadrise  angles,  the  method 
for  constructing  the  free  surface  shape  will  have  to  be  modified.  It  is 
obvious  that  the  wave  rise  cannot  be  estimated  correctly  by  Equation  (3). 
Figure  XI  is  a  sketch  of  a  typical  U-shaped  section.  The  pressure  distri- 
bution for  such  a  section  is  satisfactorily  predicted  by  the  expanding  plate 
theory,  except  in  the  spray  root  region,  from  y  »  b  to  y  ■  c.  ft  is  there- 
fore reasonable  to  expect  that  Hie  pressure  on  the  flat  portion  of  the  hull 
can  be  given  by  assuming  the  pressure  distribution  for  a  flat  plate  from 
y  ■  0  to  y  •  b.  The  flow  past  this    flat  plate    can  be  used  to  estimate  the 
wave  height  at  the  free  surface,  From  ( 18)  the  distribution  of  the  normal 
velocity  for  a  flat  plate  of  width  2a,  at  points  y  >  a*  is  given  by 
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VQ  will  b«*  the  velocity  of  the  trm  surface  particles.  The  wave  height 

at  A  poini  y  can  be  found  by; 

-t 


For  the  wave  height  determinations,  the  dimension  a  la  taken  to  he 
equal  to  cv  the  width  of  the  section  immersed  at  the  original  water 
line  (see  Fig.  XI).  When  constructing  the  free  surface,  the  spray  is 
assumed  to  shoot  off  at  an  angle  determined  by  the  slope  of  the  hull  in 
the  spray  root  area.  This  is  the  same  as  assuming  that  the  side  of  the 
ship  is  straight,  and  that  th^t  impacting  body  is  roughly  equivalent  to  a 
blunt  nosed  wedge.  Assuming  that  the  similarity  principle  still  holds, 
the  free  surface  is  constructed  by  using  the  arc  length  and  continuity 
conditions,  and  checking  by  Equation  <•). 

For  a  U-shaped  section,  the  velocity  potential  for  the  flat  portion 
of  the  hull  {from  w  •  0 to  w  »  h)  em  be  assumed  to  be  that  tor  irrotational 
flow  about  a  flat  plate.  From  (IS)  the  complex  velocity  potential  past  a 
plate  perpendicular  to  the  flow  is  given  by; 


*   -  a.  v/a*-  y8  (88) 

■ 

Equation  (if)  can  be  used  to  compute  the  velocity  potential  on  the  hull 
from  y  «  b  to  y  »  c.  The  two  velocity  potentials  should  agree  at  y  *  b 
for  any  specific  time.  (A  sample  calculation  is  discussed  in  Appendix  A. ) 

The  calculations,  from  which  the  results  in  this  paper  were  de- 
termined* were  carried  out  for  the  Lewis  type  ship  discussed  previously 
(see  Figure  XII). 


31  - 


Idglari  »vwr  eilT    »»l*l**ae;  aaaHaa  aenl  «^  to  ^liaolav  ad*  *d  lilw  flV 

:^d  fcawo*  »d  nao  \.  laioq  a  la 


ad  ol  tv*jimt  ai  *  ncliAPinlb  ad*  >aj*lfraataaiia*aa)  tdgiad  avaw  ad*  10I 

?****-  i«cJf^H<..-  cMit  I*  kaaia.rani  aoi***f)  Mil  )o  dlblw  »rft  «*9  ©I  lavpa 

al  t*taja  *&  jto&tu*  aail  nill  gaitoiralaooo  nadW  .(IX  .jl*  aaa)  •oil 

ai  JUad  <Mfl  to  aqole  »4t  t*  fcaaJnrtaiak  alfcoa  mHVd  tooda  ol  aaawea* 

•A  to  abta  oxil  tadi  galxraaaa  M  a«aa  ad*  «i  aldT  .aru  loo*  ynq*  ttft 

*  al  tn&i*v/i/f>»  'iliiguQar  ai  ^©d  saisoaqaii  adf  JtMll  baa  ,t%fmi  aj  qida 

,afefod  lilts  alqloateq  ^rtnlimi*  »dl  lull  ^alaiuaaA  ,*febew  kMi  laiild 

^iufiiJnoT*  baa  dJfeaai  on*  adf  ga/aa  ^gd  betaa -llama  al  aaaHua  —>A  *dt 

.{•)  noliaupS  td  saidaado  Mat  ,cftoJttlbfloa 

oolfvoaj  tail  aril  «sol  Uitoaleq  ^liaoiav  eatt  .aallaaa  h»aadi-U  a  «a»*V 
laaoilatoiiu  no!  iadl  do  of  MMsuaa*  9d  aao  (d  •  x  o*  •  *  t  awwt)  UuA  adl  le 
a  laaq  laiaoaloq  v;i*oolav  awlqaaao  adf  W)  mm%  .»iaiq  *aft  a  Ivoda  well 

^d  rarrtt  **  woft  adt  at  naiaaaViaejiaa;  alalq 


V  Jr  -    -     *a  -  *tVv*  "    •    * 


liiiri  ad*  iso  I«i*«ayoq  tfiooiev  ad*  etaajaaii  o*  feaav  ad  oaa  (TI ) 

d  *  \:  fa  aanja  itiaada  afalfaafoq  xnoolav  owl  adT  .a  •  \  ol  d  •  x  aiwrt 

( .A  wiwanajA  ai  feeaaitaaifc  ai  oelaaiaatea  ala/nai  A)  .ami!  oittaaqa  xp*  "?©* 

-a*  ataw  *ao>aq  aid!  ol  atfaan  adi  daiitar  mot*  .aooitalaalao  adT 
%Itu#ohraiq  baaauna  tfo  qida  awt*  olwaU  adl  ?ol  lite  oalrtaa  anew  Jfraattnwa* 

AVX  aiua**  aaa) 


-  It- 


O     I        ^        3    4         5      (,    1  11 


kO    iq       11        n         /£    IS*  IHr  H IX 


FIGURE  XII   LEWIS  SECTION  SHAPES 
-32- 


III.    RESULTS 


Included  in  this  section  are  the  curves  of  the  theoretical  pressure 
distribution  for  the  ©hip  inder  consideration.  Thees  curves  Are  shown 
In  H  gores  XIII  ■  XV  ,  Also,  one  transverse  section  of  the  Liberty  ship, 
used  as  a  model  in  reference  (3),  was  investigated.  The  pressure  distri- 
bution for  this  section  is  shown  in  Figure  XVII.  This  section  was  quite 
V-shaped  and  had  a  deadrise  angle  of  $    30*  .  Straight  wedge  theory  was 
used  for  calculation  of  the  pressure  distribution,  and  the  results  are  com 
pared  with  the  results  given  in  reference  (3)  for  the  same  section. 

A  sample  calculation  for  station  3  of  the  Lewis  section  ship  is  dis- 
cussed in  Appendix  A.  The  tangential  velocity  (v  >  on  the  flat  portion  of 
the  hull  was  calculated  by  twc  methods  (Table  IV),.  The  first  method  was 
by  the  assumption  of  flow  over  a  flat  plate.  This  is  the  method  used  for 
all  of  the  low  deadrise  sections  under  consideration.  However,  for  com- 
parison, the  tangential  velocity  on  the  bottom  of  the  hull  was  calculated 
by  use  of  the  Cauchy  integral  formula  and  equation  (17).  The  results 
check  closely  with  those  from  the  flat  plate  flow  assumption. 

The  pressures  load  on  the  forefoot  was  determined  by  strip  theory 
for  several  instants  of  time.  The  integrated  load  distribution  is  shown 
in  Figure  XVI. 
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IV    DISCUSSION  OF  RESULTS 

A,  Introduction 

Before  discussing  the  results,  the  theoretical  assumptions  made 
will  be  restated. 

1)  Only  two  dimensional  sections  were  considered. 

2)  Each  impacting  section  was  assumed  to  be  perfectly  rigid,  and  to 
be  symmetrical  about  the  center  line  plane. 

3)  While  the  velocity  of  penetration  was  not  constant,  for  the  small  time 
interval  under  consideration,  it  was  considered  to  vary  linearly  and  aver- 
age values  were  used. 

4)  The  flow  field  is  assumed  to  be  lrrotational. 

5)  The  fluid  is  considered  to  be  frictioniess  and  incompressible,  and 
the  effects  of  gravity  and  surface  tension  are  neglected. 

6)  On  the  basis  of  irrotational  flow,  the  particles  of  the  origins!  surface 
are  considered  to  remain  in  the  free  surface.  The  surface  deflected  by 
the  impacting  body  is  assumed  to  have  a  constant  arc  length. 

7)  The  theoretical  procedure  was  based  on  a  mathematical  model  of  a 
wedge  with  constant  entry  velocity,,  The  theory  was  applied  to  the  water 
entry  of  actual  ship  sections  by  assuming  that  the  impacting  sections 
were  approximately  blunt  nose  wedges.  That  is,  the  free  surface  formed 
during  the  initial  impact  stages  was  considered  to  be  equivalent  to  that 
formed  by  a  wedge  of  an  equivalent  deadrise  angle.  The  spray  was  con- 
sidered to  shoot  off  at  an  angle  determined  by  the  slope  of  the  hull  at  the 
wetted  width.  For  the  hull  sections  where  the  deadrise  angle  was  approxi- 
mately »ero,  the  velocity  potential  was  determined  by  the  expression  for 
flow  along  a  flat  plate, 
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8)  The  ship  motion  expressions  were  determined  experimentally  for  a 
Liberty  ship  of  approximately  the  same  dimensions  as  the  ship  found 
from  Lewis  sections „ 

B.  Determination  of  Pressure  Distributions 

The  ship  motion  expression  for  the  relative  velocity  did  not  include 
the  effect  of  ship  forward  speed.  For  a  small  pitch  angle,  the  total  re- 
lative normal  velocity  of  the  keel  is? 

V    ~VY    +    xi  +  *  *  r 
n 

where  U      is  the  effect  of  the  forward  speed  on  the  normal  velocity. 

If  the  pitch  angle  is  small  at  the  instant  of  slamming,  0      can  be  safely 

neglected.  At  t  *  1.284  sec,  the  first  time  for  calculations,  the  pitch 

angle  was  approximately  0.2*.  Therefore,  the  contribution  of  U      to  the 

relative  velocity  was  not  important. 

Figures  W.I  to  X  V  represent  the  pressure  distribution  on  the  fore- 
foot of  the  ship  at  approximately  the  time  of  slamming.  The  pressures 
were  computed  for  all  stations  from  1  to  5  which  are  immersed  at  the 
particular  time.  The  pressure  distributions  are  not  really  sufficient 
to  give  a  definite  value  to  the  maximum  pressure  to  which  the  ship  is 
subjected. 

The  pressure  distribution  carves  are  similar  to  the  curves  ob- 
tained in  reference  (3)  for  the  Liberty  ship.  While  the  peaks  on  the 
carves  obtained  in  this  paper  were  lower,  this  is  probably  not  significant. 
The  peak  pressures  are  of  extremely  short  duration,  and  are  relatively 
unimportant.  For  example,  at  t  *  1,270  station  3  was  out  of  the  water, 
but  .0184  sec.  later  it  had  a  peak  pressure  of  530  psi  at  y  *  5.  8?  ,  and 
.016  sec  later  there  was  a  pressure  of  approximately  120  psi  at  y  *  5.  8* 
However,  it  does  seem  that  the  average  pressure  distribution  at  a  station 
varies  directly  as  the  peak  pressure. 
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An  attempt  was  made  to  determine  the  pressure  at  station  3  at 
a  slightly  earlier  time,  but  the  construction  of  the  free  surface  was 
extremely  difficult  to  perform,  and  the  author  had  no  reason  to  believe 
the  results  were  correct.  It  is  quite  possible,  however,  that  the  peak 
pressure  might  be  even  higher  at  station  3  at  an  earlier  instant  of  time. 

Referring  to  Figure  XVII  the  pressures  on  station  %  of  the  ^Liberty 
ship  were  uniformly  lower  than  those  determined  In  reference  (3)  by 
the  expanding  plate  theory.  This  result  is  probably  the  more  accurate, 
since  the  ship  section  could  accurately  be  represented  by  a  wedge  during 
the  time  of  initial  entry  into  the  water.  The  wedge  assumption  would  not 
have  been  valid  at  a  draft  of  6  ft. ,  or  greater,  but  by  this  time  the  im- 
pact pressures  have  been  greatly  reduced,  From  the  results  given  in 
reference  <7)a  it  is  to  be  expected  that  the  computation  of  pressure 
distributions  on  wedge-shaped  bodies,  by  considering  the  free  surface, 
will  give  somewhat  lower  pressures  than  those  calculated  by  the  ex- 
panding plate  theory. 

The  integrated  slamming  load  distribution  in  tons/ ft,  of  length, 
is  shown  in  Figure  XVL  At  best,  these  curves  are  crude  approximations, 
since  each  curve  was  determined  by  only  3  or  4  points.  A  much  better 
method  would  be  to  compute  the  pressure  distribution  at   -5    and  -?  sta- 
tions in  the  area  of  slamming.  This  was  not  done  in  this  paper  since 
the  offsets  for  fractional  stations  were  not  available.  However,  even 
with  the  points  available,  a  good  idea  can  be  obtained  as  to  the  total 
slamming  load  on  the  forefoot  of  the  ship  at  the  indicated  times. 
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V.    CONCLUSIONS 

1)  The  method  of  analysis  given  in  this  paper  is  at  least  a  reasonable 
approximation  since  the  curves  of  pressure  distribution  are  of  the  cor- 
rect shape. 

2)  The  procedure  for  determining  the  pressure  distribution  by  con- 
sidering the  free  surface  shape  is  a  long  and  tedious  process.  There 
are  inherent  sources  of  error  in  the  construction  of  the  free  surface, 
and  in  the  graphical  and  numerical  solutions  for  the  pressures,  even 
for  a  simple  wedge.  The  procedure,,  when  adapted  to  a  ship  form,  is 
even  more  questionable.  However,  there  does  seem  to  be  some  merit 
in  the  procedure,  and  it  is  concluded  that  it  can  be  refined  considerably. 
The  biggest  improvement  would  be  to  develop   a  computer  solution  for 
the  procedure, 

3)  The  wave  rise  is  probably  accurately  given  for  an  impacting  body 
by  either  equation  (3)  or  (37),  depending  on  the  shape  of  the  body.  How- 
ever, it  is  evident  that  equation  (8)  cannot  be  used  to  determine  the  spray 
root  thickness  for  a  ship  section.  When  the  free  surface  was  constructed 
for  the  blunt  nosed  sections,  the  spray  root  thickness  had  to  be  estimated. 
This  was  possibly  a  large  source  of  error. 

4)  The  determination  of  the  free  surface  shape  in  the  spray  root  area 
is  quite  critical,  This  is  where  the  particle  velocities  are  the  highest, 
and  even  a  small  error  in  the  shape  results  in  a  large  mrror  in  directions 
and  magnitude  of  particle  velocities.  However,,  eveo  with  the  possible 
sources  of  error  in  the  spray  root  region,  the  pressures  actually  obtained 
for  this  region  were  as  realistic  as  those  determined  by  expanding  plate 
theory. 
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5)        The  procedure  described  in  this  paper  is  only  valid  during  the 
initial  impact,  and  a  short  time  afterward.  After  the  body  has  peae- 
trated  to  a  sizeable  depth,  the  effect  of  gravity,  and  spray  dissipation, 
will  render  consideration  of  the  free  surface  useless.  However,  the 
highest  pressures  will,  of  course,  occur  at  the  time  of  impact  or  imme- 
diately afterward,  when  the  added  mass  of  the  body  is  changing  most 
rapidly. 
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2)        Experimental  and  additional  theoretical  work  should  be 
to  determine  if  similarity  principles  can  be  applied  to  ship  sections  of 
general  shape,  Particular  attention  should  be  given  to  determining  how 
the  expansion  of  the  free  surface  is  determined  by  the  ahape  of  the  im- 
pacting body,  As  a  start,  experiments  could  be  devised  which  utilized 
impacting  wedges  of  various  deadrise  angles  and  a  high  speed  camera. 
The  spray  root  thickness  and  the  angle  at  which  the  spray  leaves  the 
body  should  be  thoroughly  investigated. 

2)  Computer  procedures  should  be  developed  for  the  solution  of 
hydrodynamic  impact  problems  involving  a  noa-ssero  velocity  potential 
free  surface, 

3)  Accurate  theoretical  and  experimental  methods  should  be  developed 
which  will  give  the  total  slam  load  acting  on  the  forefoot  of  the  ship  at 
any  instant  of  time, 

4)  The  effects  on  slamming  loads  of  hull  flexibility  and  water  com  - 
pressibUity  should  be  investigated  both  theoretically  and  experimentally. 
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APPENDIX    A 
Sample  Calculations  U -Shaped  Section  (Lawla  Ship,  Station  3), 

t  •  I.  Hi  tn, » 

The  first  estimate  of  the  surface  shape  is  computed  by  using 
equations  (37),  and  the  continuity  principle.  From  the  plotting  of  the 
wave  rise,  the  wetted  width  (c)  is  determined  to  be  5.85  ft.  (See  Figure 
A~Z),  The  slope  of  the  hull  section  at  y  *  c  is  determined  to  be  approxi- 
mately SS*  .  This  is  the  angle  at  which  the  spray  is  assumed  to  shoot 
off.  The  spray  tip,  0  •'  ,  is  located  by  measuring  along  the  wave  and 
wedge  from  y  *  10'  a  distance  a  *  y  «  101  .  (At  the  distance  y  *  15'  , 
the  wave  slope  is  quite  small,  so  that  s  »  y  is  a  close  approximation. ) 
A  spray  thickness  of  £fc  *»  At  was  assumed.  Once  the  spray  and  wave 
have  been  merged  together,  the  spray  tip  is  located  more  exactly  by 
following  the  curbed  surface  ($'  is  located  or*  the  line  drawn  tangent 
to  the  hull  at  y  s  c,  at  an  angle  of  f!  •  S3*).  The  area  of  the  displaced 
water  above  the  original  water  line  is  made  to  equal  the  area  of  the 
submerged  body  below  the  original  water  line.  Then  equations  (6)  are 
used  to  check  the  position  of  the  water  particles,  The  steps  of  the  num- 
erical integration  for  the  surface  shape  are  shown  in  Table  II,  The 
new  particle  positions  were  determined  by  the  following  approximate 
forms  of  equations  (6)i 
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After  the  free  surface  shape  has  been  determined,  the  velocity 
potential  can  be  determined  on  the  body.  For  the  flat  portion  of  the 
hull,  from  y  *  0  to  y  *  4. 7,  the  potential  is  determined  by  equation  (28). 
In  this  computation,  a  is  taken  to  be  equal  to  the  wetted  width  at  the 
time  under  consideration,  but  the  equation  is  only  used  for  the  flat 
portion  of  the  hull  (from  y  «  0  to  y  *  4.  7). 

For  the  portion  of  the  hull  {from  y  *  4. 7  to  y  *  6*.  96)  the  velocity 
potential  is  computed  by  integrating  the  velocity  along  the  hull  surface. 
The  tangential  velocity  (v  >  along  the  hull  is  determined  from  equation 
(17).  Four  points  are  chosen  on  the  curved  portion  of  the  hull  for  the 
determination  of  the  velocity  (see  Figure  A«i).  The  iy  and  x  axis  are 
constructed  tangent  and  perpendicular  to  the  buiU  at  each  point  under 
consideration.  The  velocity  <v  >  is  determined  by  equation  ( tf )  by  num- 
erical integration  around  the  free  surface; 

v    *  — \   vd&  +  -  S    u  d  (la*) 

The  steps  in  the  integration  are  shown  in  Table  III  for  point  2  (j,,  vf  and  u 
refer  to  the  position  snd  velocities  of  points  on  the  free  surface),,  On  the 
body,  v    is  equal  to  s    cos  3  at  any  point.  Knowing  the  component  velo- 
cities,  the  resultant  velocity  (U)  can  be  found  at  any  point.  The  velo- 
cities and  velocity  potential  determined  at  the  points  on  the  hull,  shown 
in  Figure  A -I,  are  listed  in  Table  IV. 

Now  the  pressure  can  be  determined  for  equation  (34); 

p        d>       t  o;- 

-E_         — 1_    4-   —    (j  £os  a  *•    -r- 
P      •      1  t  2 

where     is  the  distance  from  the  origin  to  the  point  on  the  body  under 
consideration,  and  a  is  the  magi®  between    and  0.  The  results  of  DM 
solution  of  equation  (24)  for  station  •  at  t  •  i.284  are  listed  in  Table  IV. 
The  time,  t,  involved  in  equation  (24)  is  the  elapsed  time  from  impact 
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of  the  body.  For  the  caae  under  consideration,  t    (impact  time)  ■  1.270, 
and  the  computation*  were  carried  out  from  time  t*    »  1. 284  sec.  There - 
fore,  I  •%•  I    "  .014  aec*  in  equation  (34). 
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TABLE  II 


Station 


NUMERICAL  INTEGRATION  FOR  SURFACE  SHAPE 
3  t  =  1.284  zx  =  0.3  zx   =  26.0  $   =  63° 


s 

V 

wt 

wt 

vrt  A  s 

7? 

• 

start 

s3- 

s* 

new 

o< 

Ut 

10 

.13 

.15 

.001     I 

.0025 

.013 

.13 

90 

.15 

9 

.14 

.30 

.004 

.007 

.016 

.14 

90 

.30 

8 

.17 

.60 

.009 

.015 

.022 

.18 

83 

.60 

7 

.24 

1.00 

.021 

.017 

.037 

.25 

85 

1.00 

6.5 

.31 

1.90 

.045 

.013 

.054 

.35 

79 

1.90 

6.25 

.45 

4.00 

.103 

.026 

.080 

.50 

69 

4.20 

6.05 

.59 

5.80 

.154 

.026 

.106 

.64 

55 

6.70 

5.9 

.71 

6.50 

.187 

.037 

.132 

.77 

45 

8.30 

5.7 

.89 

6.00 

.134 

.038 

.169 

.96 

27.5 

10.55 

5.5 

1.09 

6.00 

.198 

.207 

1.15 

25.3 

10.55 

Area  of  submerged  body  =  1.54  ft 
Area  of  displaced  water  =  1.60  ft 


10 
9 
8 
7 

6.5 
6.25 
6.05 
5.9 
5.7 
5.5 


start 


0 

0 

0 

.04 
.07 
.09 
.13 
.30 
.53 
.80 


vt 


0 
0 
0 
.08 

.25 
1.30 
3.30 
5.20 
8.70 
8.70 


vt 


0 

0 

0 

.0016 

.006 

.033 

.090 

.149 

.268 

.288 


vtAs 

S** 


•vtz^s 


.004 

.005 
.012 
.018 
.042 
.056 


0 

0 

0 

0 

.004 
.009 
.021 
.039 
.081 
.137 


Tic 

new 


0 
0 
0 
0 

.03 

.06 
.13 

.24 
.46 
.75 
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TABLE  III 
Numerical  Integration  of  Equation  (17) 


z 

V 

u 

e 

vdG 

R 

In  R 

ud(ln  R) 

s  =  10 

10 

-  5 

.40 

0 

4.2 

1.43 

8.0 

s  =  8 

38 

-  20 

.40 

-  8. 

2.2 

.79 

33.0 

s  =  6.5 

129 

-  42 

.30 

163 

.75 

-   .29 

3.4 

s  =  6.05 

475 

25 

.84 

245 

.50 

-  .69 

34.0 

s  =  5.70 

666 

354 

1.27 

93 

.60 

-  .51 

124 

s  =  5.50 

666 

354 

1.41 

80 

.85 

-  .16 

478 

s  =  3.0 

666 

354 

1.53 

27 

3.3 

1.19 

230 

s  =  0 

666 

354 

1.57 

600 

910 

v0  (point  2)  =  ^(600)  +  7f(910)  =481  ft/sec 
Uo  =  zx  cos63  =  11.8  ft/sec 


TABLE  IV 

Velocity,  Velocity  Potential,  and  Pressure  Distribution 
on  Station  3  of  Lewis  Ship(t  =  1.284  sec) 


y 

0 

v0 

Vo 

u0 

P 

-  . . 

(Plat  Plate) 

(Cauchy  Integral) 

(f*0 

0 

-154.5 

0 

5 

26.0 

144 

1 

-152.0 

4o43 

6 

26.0 

145 

2 

-143.5 

9.42 

10 

26.0 

147 

3 

-134.0 

15.16 

17 

26.0 

165 

4 

-114.8 

23.60 

24 

26.0 

200 

4.7 

-  98.2 

32-30 

28 

26.0 

202 

5.6 

3.0 

• 

158 

20.5 

550 

5.85 

160 

481 

11.8 

200 
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